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1. THE PROBLEM OF THE DETERMINED ANTS. Imagine four determined 
ants who simultaneously walk along the edges of the picnic table graph of Figure 1. 
The ants can move only to the right (northeast, southeast, and sometimes due east) 
with the goal of reaching four different morsels. 




Figure 1. Ants and morsels. 



Question 1. How many ways can the ants in Figure 1 simultaneously reach different 
morsels? 

See Figure 3 for an example. We define an n-path to be a collection of n paths from 
a set of n origins to a set of n destinations. To compute the number of 4-paths in our 
example, we first find the number of ways that each ant can reach each morsel. In the 
example these numbers can be computed easily using calculations similar to those that 
arise in Pascal's triangle (see Figure 2). We record the information in a square matrix 
A whose {i, 7) -entry aij is the number of ways that Ant ; can reach Morsel j. 



Ant 2 




1 4 

Figure 2. The number of paths for Ant 2 to reach each morsel can be computed recursively and is recorded in 
the second row of matrix A. 
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Consequently, the number of 4-paths such that each Ant / reaches Morsel / is 
an'322'333a44 = 44100. Alternatively, for each permutation it of {1, 2, 3,4} there are 
'3i3r(i)'323r(2)'333r(3)<24jr(4) 4-paths whcrc cach Ant / reaches Morsel 7T{i). Summing over 
all alternatives gives us 

Ol;r(l)<22;r(2)03jr(3)fl4;r(4) 

4-paths, where ^4 is the set of all permutations of {1, 2, 3, 4}. In other words, the 
"antswer" to Question 1 is the permanent of the matrix A. In our example, this equals 
171361. 

Question 2. How many ways can the ants simultaneously reach different morsels, 
where no two paths intersect? 

Notice that a feature of our picnic table graph is that in order to have four non- 
intersecting paths, we must have each Ant i go to Morsel /. Such a graph is called 
nonpermutable. However, most of the 44100 4-paths associated with the identity per- 
mutation do intersect somewhere. Believe it or not, the answer to Question 2 is the 
determinant of A, in this case 889. 

To see why this is true, first recall that for an n-hy-n matrix A the determinant of A 
equals 

sgn(7r)ai;r(i)'223r(2) ■ ■ ■ a„„(„), 

1T€S„ 

where 5',, is the set of permutations of {1, 2, . . . , n} and sgn(:T) = 1 when it is an 
even permutation (expressible as the product of an even number of transpositions) and 
sgn(7r) = — 1 when tt is odd. In terms of our counting problem, the determinant is 
a weighted sum, over all n-paths, where the weight of an «-path from Ci, . . . , C„ to 
respective destinations D^^d), . . . , D^f,,-) is the sign of it. In a nonpermutable graph, 
all of the nonintersecting «-paths are associated with the identity permutation, which 
is even, and thus given positive weight. It remains to prove that for every intersecting 
n-path, we can uniquely identify another intersecting n-path with opposite sign. 

For a given intersecting n-path with associated permutation jt, suppose / is the 
smallest index for which Path i intersects another path, and let Path j be the largest 
indexed path that Path / intersects. Let O be the first point of intersection of Paths i 
and j. For the 4-path in Figure 3, ; = 1 and j = 4. To create an intersecting n-path 

.4 .4 




Figure 3. An intersecting 4-path with even permutation n = (13) (24) is transformed into another intersecting 
4-path with odd permutation n' = (1243). 
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with opposite sign, we simply swap Paths / and j after O. Thus, Path / ends up at D,^^j), 
Path j ends up at Z)jr(,), and all other paths remain unchanged. Hence, the new n-path 
will have associated permutation n' = (i, j)jr, which necessarily has opposite sign. In 
Figure 3, tt = (13)(24) is even and n' = (14)(13)(24) = (1243) is odd. Notice that 
in the new n-path j, and O are the same as before, so (tt')' = tt, and the process is 
completely reversible. 

The preceding argument applies to any directed graph that is acyclic (i.e., has no 
cycles) and leads to the following theorem: 

Theorem 1. Let G be a directed acyclic graph with n designated origin and destina- 
tion nodes, and let A be the n-by-n matrix whose (/, j)-entry is the number of paths 
from the ith origin to the jth destination. The following statements hold: 

(a) The number of n-paths is equal to the permanent of A. 

(b) If G is nonpermutable, the number of nonintersecting n-paths is equal to the 
determinant of A. 

(c) In general, even if G is not nonpermutable, the determinant of A equals 
Even{G) — Odd(G), where Even(G) is the number of nonintersecting n-paths 
corresponding to even permutations and Odd(G) is the number of noninter- 
secting n-paths corresponding to odd permutations. 

This theorem was originally given by Karlin and McGregor [5] and Lindstrom [7], 
and popularized by Gessel and Viennot [4], Aigner [1], and the recent book by Bres- 
soud [2]. When the given graph G has considerable structure (as in cases where it is 
used to enumerate Young tableaux, plane partitions, or rhombus tilings), it is often 
possible to find a closed form for the determinant of A. Krattenthaler has described 
several methods for evaluating such determinants in [6]. 

Application. Recall that for n >Q the nth Catalan number 

1 /2n\ 
n-\- \\n J 

counts the lattice paths from (0,0) to {n,n) restricted to vertices that stay on or 
below the line y = x. Let H be the {n + l)-by-(n + 1) matrix with (;', 7')-entry 
Ci+j (0 < (, j < n). In Figure 4, the number of paths from origin O, to destina- 
tion Dj is Ci+j. Since there is only one way to create 7? + 1 nonintersecting paths 

D2 



Oo = Do 



O2 



O3 



Figure 4. There is only one nonintersecting 4-path from {Oq, Oi, O2, O3) to {Dq, D\, D2, fj). 



det 



" Co Ci C2 C3 " 

C\ C2 C3 C4 
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from {Oo, Oi, . . . , O,,} to {Dq, Di, . . . , D,,}, Theorem 1 implies thatdet(//) = 1. For 
another quick combinatorial derivation using matrix factorization, see [12]. 

2. ENUMERATING SPANNING TREES. Now consider the graph G in Figure 5. 
A spanning tree of G is a connected acyclic subgraph of G with the same vertex set. 



10 10 




3 7 8 3 7 8 

Figure 5. A graph G (left) and one of its many spanning trees (right). 



Question 3. How many spanning trees does G have? 

To answer this question, we first describe the graph in matrix notation. We allow 
our graph to contain multiple edges (e.g., between vertices 7 and 8 in G) but no loops 
(edges that begin and end at the same vertex). For a graph with n vertices, let A be 
its n-hy-n adjacency matrix, where a,-,- is the number of edges between / and j. Let 
D be the diagonal matrix Diag((ii ,d2, . . . ,d„), where J, is the degree of vertex / (i.e., 
the number of edges incident at i). Now consider the matrix D — A. Since every row 
of D — A sums to zero, its columns are linearly dependent; hence, det(D — A) = 0. 
However, if we remove any row and column from D — A, the determinant of this 
submatrix answers our question. Specifically, we have: 

Theorem 2. Let G be a loopless undirected graph with n vertices, adjacency ma- 
trix A, and diagonal degree matrix D. If B,s signifies the {n — V)-by-(n — 1) matrix 
obtained by deleting from D — A its rth row and sth column, then the number of span- 
ning trees of G is equal to (— 1)' det B, s for any choice of r and s. 

This result is known as the "Matrix-Tree Theorem" and was originally proved by 
Sylvester [15] (see [17], [14] for algebraic proofs). We present a combinatorial proof 
(as given in [3]) of the case when r = s = n, whereupon the number of spanning trees 
is simply the determinant of B,,,,, a matrix that we abbreviate as B. 

For example, the graph in Figure 5 has 
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The 9-by-9 matrix B obtained by removing the last row and column from D — A is 
seen in Figure 10 (with its nonzero entries expanded). It has determinant 148. Hence 
the matrix-tree theorem tells us that G has 148 spanning trees. 

Our proof of Theorem 2 is similar to the one in the previous section. Here, we will 
identify spanning trees as the acyclic objects in a large collection of directed graphs. 
The determinant will count all the acyclic objects and half of the cyclic objects posi- 
tively, while the other half of the cyclic objects will be counted negatively. When the 
dust settles, only the spanning trees will remain standing. 

For any spanning tree of G with n vertices, there is exactly one way to orient its 
edges so that each edge points in the direction of vertex n (see Figure 6). This is called 
a rooted spanning tree with root n. Hence, the number of spanning trees of G is equal 
to the number of rooted spanning trees of G with root n. 

10 




3 7 8 

Figure 6. The spanning tree of Figure 5 rooted at vertex n = 10. 

Observe that the rooted spanning trees of G have n — I edges, where each nonroot 
vertex / has outdegree 1. A directed subgraph of G with this property is called afunc- 
tional digraph in G, since it represents a function / : {1, — 1}^ {1, ...,«}, 
where the edge from ; to j indicates that /(/) = j. Let !F denote the set of all func- 
tional digraphs in G. Notice that the size of is d\d2 ■ ■ ■ d„-i (see Figure 7 for a typical 
example). From any vertex, we are ultimately led either to the root n or to a cycle that 
does not contain n. The acyclic objects of are precisely the rooted spanning trees 
of G. 

10 




Figure 7. A typical functional digraph F in G. 

Next, we consider an even larger collection S of signed objects: S comprises all 
functional digraphs F in G, where now each cycle in F is given a sign, either + 
or — . Thus a functional digraph F in J- with k cycles generates 2* signed functional 
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Figure 8. S is a signed functional digraph based on F from Figure 7. Here, the sign of S is negative. 

digraphs in S (see Figure 8). For 5 in 5 we define the sign sgn(5') of 5* to be the 
product of the signs of the cycles of S. If S has no cycles, then 5' is a rooted spanning 
tree and its sign is necessarily positive. If S is cyclic (that is, S has at least one cycle), 
then we define its conjugate S to be the same functional digraph as S but with the 
sign of the first cycle reversed. We define the first cycle to be the cycle containing 
the vertex with the smallest label (see Figure 9). Notice that sgn(5') = — sgn(5') and 

S = S. This correspondence proves that there are as many positive cyclic elements as 
negative cyclic elements in S. Hence the number of positive elements of S minus the 
number of negative elements gives the number of acyclic elements of S, which is the 
number of spanning trees of G. 




Figure 9. 5, the conjugate of S from Figure 8, has positive sign. 



So how does the matrix B fit into this? Suppose that an edge between / and j 
exists in G, where / ^ n and j ^ n. Then the directed edge from / to j is represented 
twice in B: positively, on the diagonal, as one of the Is that comprise ba = J, = 
1 + 1 + • • • + 1, and negatively, off the diagonal, as one of the —Is that comprise 
bij = —I — I ■ ■ ■ — I. (Note that bij = — 1 unless there are multiple edges between 
vertices / and j.) A directed edge from / to n is represented only once, positively, as 
a 1 in bjj. Thus, every signed functional digraph S, which consists of 7? — 1 directed 
edges, is associated with a selection of n — 1 Is or —Is in the matrix B, where each 
1 or — 1 comes from a different row and column of B. Specifically, for an edge eij 
in S directed from / to j that belongs to a negative cycle of S we select one of the 
— Is in bij = —I — I — ■■■ — I. Otherwise, we utilize the appropriate 1 in = J, = 
1 + 1 + ■ ■ ■ + 1 (see Figure 10 for an example). 

Thus, when we expand each bij as the sum of Is or —Is, the determinant of B is 
represented as the sum of the products of a bunch of Is, —Is, and Os. Each nonzero 
term corresponds to a signed functional digraph 5' in G. If S has exactly m edges in 
negative cycles, then its contribution to det(fi) is Xs = (—1)'" sgn(7r5), where its is 
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the permutation associated with S. Our goal is to show that Xs = sgn(5'), from which 
it follows that the determinant of B equals the number of positive elements of S minus 
the number of negative elements of S (i.e., the number of spanning trees of G). 
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Figure 10. The signed functional digraph S in Figure 8 corresponds to the selection of 9 bold Is and — Is in 
the matrix B. 



To prove that Xs = sgn(S), suppose that S contains k negative cycles Ci, . . . , Q, 
where C, contains m, (> 2) edges. Thus, sgn(5) = (—1)*^, mi + ■ ■ ■ + m^. = m, 
and Its = Tt\iT2 ■ ■ - Ttk, where tt, is the cyclic permutation with m, elements natu- 
rally described by C,. For example, in Figure 8, its = (37), whereas in Figure 9, 
TTj = (296)(37). Thus, since sgn(7r,) = (—1)'"'"', we have 

Xs = (-1)'" sgn(;rs) = (-1)"' sgn(7ri • • • tt,) 
= (-l)""+-+"'*sgn(7ri)...sgn(7r,) 
= (^—iy"i+ -+"'k(^_iy"i-i . . . (_]^y'u-i 

= (-1)* = sgn(5), 

as was to be shown. ■ 



In the foregoing proof of Theorem 2, we assumed that (= B„„) was obtained 
from D — A by deleting its last row and column. If r = s ^ n, then the theorem can be 
proved in the same way by considering trees rooted at vertex s.li r ^ s, then a more 
sophisticated combinatorial proof is provided by Chaiken [3]. Other combinatorial 
arguments are also given by Orlin [11], Temperley [16], and Zeilberger [18]. 

The matrix-tree theorem can be extended in several directions, all of which can 
be proved by essentially the same combinatorial argument. If we let 5+ be the non- 
negative matrix D + A with its last row and column removed, then the size of S 
is the permanent of since each signed functional digraph is counted once. Fi- 
nally, notice what happens if we remove the last row and column of B, resulting in an 
(n — 2)-by-(n — 2) matrix C. An argument similar to the original one then reveals that 
the determinant of C counts the spanning /ore^fj (i.e., acyclic subgraphs) F of G with 
two connected components, where n — I and n are in different components of F. In 
general, we have: 

Corollary 3. Let G be a loopless, undirected graph with n vertices, adjacency matrix 
A, and diagonal degree matrix D. For a given set of vertices V = {vi, . . . , Vk} let By 
and By denote the (n — k)-by-(n — k) matrices obtained by deleting rows and columns 
V\, . . . ,Vk from D — A and D + A, respectively. 
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(a) The determinant of By counts the spanning forests F of G with k connected 
components, where the vertices V\, . . . ,Vk are in different components of F. 

(b) The permanent of By counts the signed functional digraphs of G with exactly 
k connected components, each rooted at one of the vertices V\, . . . ,Vk- 

Application. The complete graph K,, is a graph with n vertices and (j) edges whose 
vertices are pairwise adjacent. Cayley's formula asserts that K„ has n"~- spanning 
trees. To see this, observe that each vertex has degree « — 1, so 5 = — 7 is the 
associated (n — l)-by-(n — 1) matrix, where / is the identity matrix and J is the matrix 
of all ones. Since J has rank one, it has eigenvalue with multiplicity n — 2 and 
eigenvalue n — \ (with eigenvector (1, 1, . . . , 1)). Thus, B = nl — J has eigenvalues 
n and 1, where n has multiplicity n —2. Consequently, the number of spanning trees of 
Kn is det(B) = n""^. The complete bipartite graph K^ji has m + n vertices Ui , . . . , Vm, 
wi, . . . ,w„ and mn edges, one edge for each v, , wj pair. As an exercise, we invite the 
reader to show that K„ „ has m""'n"'"' spanning trees [8]. 

3. PERMUTATIONS WITH SPECIFIED DESCENTS. In this section, we count 
restricted arrangements of numbers. The arrangement 382469157 has descents occur- 
ring in positions two and six, since the second number and sixth number are immedi- 
ately followed by smaller numbers. 

Question 4. How many arrangements of { 1 , 2, . . . , 9} are possible, with the restriction 
that a descent is allowed at positions two, six, or seven but at no other positions? 

For example, 382469157 is counted among the valid arrangements, as is 123456789, 
since it has no descents. The answer to Question 4 is the multinomial coefficient 

9\/7\/3\/2\ 9! 

= 3780, 



2J\4J\IJ\2J 2!4!1!2! 

since this counts the ways to select two numbers to occupy positions one and two, four 
numbers to occupy positions three through six, one number to occupy position seven, 
and two numbers to occupy the last two positions, with each selection of numbers 
written in ascending order. In general, if 5' = {^i ,...,>?<;} is the set of positions where 
a descent is allowed in an arrangement of {1,2, . . . ,n — 1}, then the number of valid 
arrangements is 

n\ 

f{n; S) = . (1) 

si\ {S2 - si)\ ■■■{Sk- (n - SkV. 

Question 5. How many arrangements of the numbers {1,2, . . . , 9} are possible, with 
the restriction that a descent must occur at positions two, six, and seven but at no other 
positions? 



We can easily compute the answer 
9! r 9! 9! 9! 



2! 4! 1!2! 



6!1!2! 2!5!2! 2!4!3! 



+ 



9! 9! 9! 

+ + 



7!2! 6!3! 2!7! 



9! 

9! = '''' 



using the principle of inclusion-exclusion as follows. From the set of /(9; {2, 6, 7}) 
arrangements previously considered, we subtract those arrangements without descents 



488 



© THE MATHEMATICAL ASSOCIATION OF AMERICA [Monthly 112 



in positions two, six, and seven, respectively; we add back arrangements without de- 
scents in positions two and six, two and seven, and six and seven, respectively; finally, 
we subtract the single arrangement without any descents. Each of these subproblems 
can be computed using equation (1). For example, the first subtracted term, which 
counts arrangements without descents in position two, but with possible descents in 
positions six and seven, is equal to /(9; {6, 7}) = 9!/(6! 1! 2!). 

In general, let e{n; S) count the arrangements of the numbers 1 through n with 
descents occurring precisely at positions si, . . . ,Sk. Then inclusion-exclusion gives us 



e(n;5) = ^(-l)^-l^l/(n;r). 



(2) 



As we will soon see, e(n; S) can be computed with the aid of a determinant. For 
example. 



e(9; {2,6,7}) = 9! det 



■- i 1 1 1 

2! 6! 7! 9! 



1 i: 
1 

L 



In general, we have: 



Theorem 4. For a given subset S = {si, . . . , Sf,} of {1, — 1} let C be the 

(k + l)-by-{k +1) matrix defined as follows: if i > j + I, Cjj = 0; ;/ ; < 7 + 1, 
Cij = l/(Sj — s,-i)!, where sq = and s^+x = n. Then 

e{n\ S) = n\ det(C) 

counts the arrangements of {1, ...,«} with descents occurring precisely at positions 

Si,...,Sk. 



The matrix C is almost upper-triangular with Cj+i j = 1 (1 < j < A:) on its subdi- 
agonal and only Os below the subdiagonal. Here, we have 

nl det(C) = ^nl sgn(7r)ci^(i,C2^(2) • • • Ck+i^^k+i), (3) 

wen 

where O denotes the set of all permutations of { 1 ,...,/:+ 1 } that satisfy 7t{j + I) > j 
for j = I, . . . ,k. All other permutations result in a product of 0. 

To prove Theorem 4, we prove that the summands in equations (2) and (3) are 
identical. Notice that (2) involves 2'' summands. We claim that equation (3) also has 
2* summands by finding a one-to-one correspondence between subsets of {1, . . . , ^} 
and n. Indeed, if 7 is a subset of {I, ... ,k} with complement 7', then we associate 
with / the permutation jt in n that satisfies jT(j + I) > j if and only if j belongs 
to J. In other words, it is completely determined by 7' , the columns from which we 
select subdiagonal Is. 

For example, suppose that k = 20 and J = {4, 10, 11, 18}. We construct the unique 
permutation jt that chooses the subdiagonal element c^+i,j = 1 for every column j 
(i.e., jT(j + 1) = j), except when j is a member of {4, 10, 11, 18, 21}. Since 1, 2, and 
3 are not elements of 7, we must choose the subdiagonal Is in columns 1, 2, and 3, 
which are the subdiagonal entries of rows 2, 3, and 4. Thus, from column 4, which 
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corresponds to a member of 7, we may not choose the entries of rows 2,3,4, or 5, 
for C54 is on the subdiagonal. Hence cu is the only eligible nonzero term in column 4. 
Thus, we must have 

7r(l)=4, 7r(2) = l, :r(3) = 2, ;r(4) = 3. (4) 

Since 10 is the next element of 7, we are forced to choose subdiagonal Is in columns 
5, 6, 7, 8, and 9 or, equivalently, in rows 6, 7, 8, 9, and 10. Thus, we cannot select the 
entry c,,io for / = 6, 7, 8, 9, or 10, nor for / = 1, 2, 3, or 4 by (4), nor for ; = 11, since 
10 belongs to /. Accordingly, 05,10 is the only eligible choice in column 10. We infer 
that 



7r(5) = 10, 7r(6) = 5, n(7} = 6, :r(8) = 7, 7r(9) = 8, ;r(10) = 9. 

Continuing with this logic, we must also select Cn n, Ci2,i8> 21; all other selec- 
tions are subdiagonal Is. Consequently, in this example, the unique product associated 
with the columns of J is 



nl Ci4C5^ioCn_ 11^12.18^19,21 



54] (sio - 54)! (^11 - 5io)! (ii8 - ill)! (n - si^y. 
fin; {4, 10, 11, 18}). 



Here, the associated permutation is 

TT = (1, 4, 3, 2)(5, 10, 9, 8, 7, 6)(11)(12, 18, 17, 16, 15, 14, 13)(19, 21, 20), 

which has sign (—1)^'"^ = 1. 

In general, with a subset J = {ji, ■ ■ ■ , j,,,} of {l, . . . ,k] we associate TTy in D, 
where 7ry(7 + 1) = j for j in y, and 

7ry(l) = jl, 
^jijl + 1) = 72, ■ ■ ■ , TtjUm-\ + 1) = jm, 

nj{u + \)=k+\. 

Since TTy has /:+ 1 elements inm + 1 cycles, sgn(7ry) = (— 1)*"™ = (—l)*"!-'!. Letting 
Tj = {sj^, . . . , it follows that 

n\ det(C) = 5^ sgn(7ry)«!ci,^-,c^-,+i,;, • • • c^-,„_,+i,;,„c^-„+i,^+i 

nj€U 

= y (-l)*-l^l 

= 5^(-l)*-l^^l/(n; Tj), 

TjCS 

as desired. ■ 



This theorem was originally proved by MacMahon [9]. It is interesting to note that 
Theorem 4 can also be proved as an application of Theorem 1, as demonstrated by 
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Gessel and Viennot [4]. An extension to determinants with polynomial entries (q- 
binomial coefficients) is given by Stanley [13]. 



Application. There is only one permutation of {1, 
position but the last one. That is, e{n; {1,2, ... ,n 
trix C with Cij = l/ij — i + 1)! if / < y + 1 and c,j 
n\ det(C) = 1. For example. 



. ,n} with a descent at every 
1}) = 1. For the n-hy-n ma- 
: otherwise we conclude that 



det 
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More combinatorial approaches to linear algebra are presented by Zeilberger [18]. 
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Goldbach's Conjecture Implies Bertrand's Postulate 

In recent years, there has been a great deal of renewed interest in Goldbach's 
Conjecture (GC) that every even integer greater than two is the sum of two 
primes, due primarily to the publication of a novel [1] and the substantial prize 
offered by its publishers for a proof. (The prize was not claimed by the deadline, 
and right now the best result is J. R. Chen's 1966 theorem that every sufficiently 
large integer is the sum of a prime and the product of at most two primes.) 

An examination of several number theory texts, old and new, reveals that an 
interesting consequence of GC seems to have been overlooked: GC provides a 
short proof of Bertrand's Postulate, which states that for every positive integer 
greater than 1 there exists a prime p such that n < p < 2n. The proof is simple. 
Since n > \,2n is an even integer greater than two, whence 2n = p^ + p2 for 
primes p\ and p2- If both pi and p2 are less than n, then p\ + P2 < 2n. There- 
fore, at least one of the p,- is greater than or equal to n — that is, n < p, < 2n. 
If n is not a prime, strict inequality holds. On the other hand, if n is prime, then 
« + 1 is composite and 2{n + 1) = p[ + p'2, where p[ or p'2 (say, p'2) must be 
greater than n + 1 — that is, greater than n. Accordingly, n < p'2 < 2n + 2. But 
p'2 can't be 2« + 1 (or else p'^ = I) nor can it be 2n (since 2n is composite), so 
n < p'2 < 2n. 

The usual proof (see, for example [2, pp. 367-368]) requires some reasonably 
sophisticated analytic estimates and is omitted from most elementary texts. 
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submitted by Yoshihiro Tanaka, Hokkaido University, Japan 



492 



© THE MATHEMATICAL ASSOCIATION OF AMERICA [Monthly 112 



